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' Abstract 

We introduce Koenderink niters as novel tools for statistical cosmology. Amongst several 
■ promising applications, they provide a test for the Gaussianity of random fields. We focus 

^ ■ on this application and present some preliminary results from an analysis of the Cosmic 

Microwave Background (CMB). 

\o '. 

1 Introduction 

i> : 

The anisotropies in the Cosmic Microwave Background Radiation (CMBR) are the oldest features 
of the Universe accessible to observations. Coming directly from the last scattering surface, they 

Qh| supply information on structure at early epochs. Therefore they greatly help in the cosmologists' 

q ■ task to constrain and outrule various models of structure formation. 

^ ■ However, the Microwave sky as observed today is not merely of cosmic origin. Various sources 

of noise and other signals obscure and distort our observations. It is thus essential to apply 
sophisticated data analysis techniques in order to gain at least some insight into cosmic evolution. 
To overview just a few of the many methods employed, consult Kogut et al. (1996), [Tegmark 
(1996)| , or [Ferreira fc Magueijo (1996)| . 

In this talk we present a novel statistical method that originated in medical imaging. Considering 
the temperature anisotropies as a realization ^o(x) of a continuous random field, we systematically 

• remove structure on specific scales, and thus enhance signal over noise, 

• describe local geometry by a small yet complete set of measures, 

• construct meaningful and interpretable descriptors to assess statistical properties. 



Along these lines of thought, [Koenderink (1984)| developed an approach to the analysis of flat, 



two-dime nsional images that is n ow successfully implemented in computer vision and medical 
imaging ( |ber Haar Romeny 1996| ). We have generalized his approach to fields both on curved 
supports and in higher dimensio ns. A number of applications to cosmological datasets have been 
developed in Schmalzing (1996)] , we will now briefly outline the method and then present some 



preliminary results from an analysis of the CMB. 
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Figure 1: We demonstrate the idea behind the scale-space operator Q a described in the text. 
Structure on scales below a is removed from the original field i^o, leaving a new field u a . In this 
example, the separation of the two spikes disappears as filtering with a larger scale merges them 
into one featureless blob. 

2 The Method of Koenderink Filters 

2.1 Scale-Space Operators 

In order to separate the various scales present in a given field uq we wish to erase structure on 
scales below a given length a, thus arriving at a new field u a . In the most general approach, this 
is done by applying an operator Q a: that is 

U a = GcjUq. (1) 

The precise form of the so-called scale-space operator remains to be specified. 

It has been shown by Koenderink fc van Doom (1987)| (see also the appendix of ter Haar Romeny 



et al. 1991) that the natural choice for Q a is the convolution with a Gaussian kernel of width a: 

^,y) = I ^e x p(!!i^). (2) 

In order to achieve this uniqueness we require the following simple and actually fairly compelling 
properties of the operator: 

• Additivity and Closure: Iterating two scale-space operators is equivalent to applying a 
single scale-space operator with an "added" scale; it turns out that only with additivity of 
the squares, that is g\ gi = \J o\ + cf\ we can sa tisfy all requirements simultaneously. 

Sa, O g a2 = (3) 

• Limits: Erasing no structure leaves the image unaltered, while erasing structure on all scales 
erases the image completely. 

g a ^o = 1 and Ga^oo = 0. (4) 

• Linearity: We assume validity of the superposition principle, that is structure on one scale is 
not influenced by features on another scale. The operator is thus linear and can be expressed 
as a convolution with an appropriate kernel, 

w CT (x) = J d d x J d rf yG CT (x,y) Wo (y). (5) 
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Figure 2: These panels illustrate the application of scale-space operators (see Section pj] ) and 
geometric invariants (see Section |2.2|) to a two-dimensional image. All nine panels were constructed 
from the image in the top left corner. The first row shows a pentagon and two of the geometric 
invariants, namely edgeness and cornerness which obviously enhance the corresponding features. 
In the second row we see the pentagon disturbed by Poissonian noise points of equal strength. 
Although the pattern is still visible, the geometric invariants - again we show cornerness in the 
middle column and edgeness in the right column - fail to find the inherent features. However, by 
applying a scale-space filter of moderate width, we arrive at the series of images shown in the 
bottom row. Both the original pentagon and its edges and corners as seen by geometric invariants 
have been reconstructed, although the noise is still present. 
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edgeness 

profile's extrema . . . 
. . . and saddlepoints 
cornerness 



Table 1: This table summarizes a choice of basic invariants that is unique and complete to 
second order in two dimensions. Furthermore, all invariants may be interpreted as morphometric 
quantities as given in the right column. 



Invariance: Rotational and translational invariance reduce the dependence of the kernel 
on two points to a dependence on a single scalar, the distance of the points. 



G a (x,y) = G CT (||x-y||). (6) 



2.2 Geometric Invariants 



The simplest way to characterize local geometry around a point of the support of a random field is 
to expand the field as a Taylor series to sufficient order N around this point. The first two orders, 
and the general expression are^] 

?v(x + y) = u ayi (x.)yi + 0(y 2 ) 



u a (x + y) = u aii (x)yi + ^<7,ij(x)y*2/j + 0(y 3 ) 



(7) 



M=0 

Obviously, considering the derivatives up to N-th order allows us to gain insight into the field's 
local behaviour while discarding higher orders in a controlled manner. However, mere derivatives 
depend on the choice of coordinate system. If we remove these unwanted degrees of freedom by 
considering coordinate independent invariants such as the square of the gradient instead of its 
components, we arrive at a description of local geometry that is both complete, in the sense that 
it contains all information up to N-th order, and unique, in the sense that all invariants can be 
constructed as polynomials in a few basic invariants. 

As an important example we will consider a basis for invariants up to second order in two dimen- 
sions. By counting degrees of freedom - derivatives of first and second order give five independent 
numbers, but one of these degrees of freedom is removed when fixing the coordinate system - we 
conclude that four invariants, apart from the field u itself, can be found. A simple choice is 

(8) 

Not all of these basic invariants allow for an intuitive interpretation. We have to go to a different, 
equivalent basis which is summarized in Table together with geometric interpretations, u^u^ is 
the square of the gradient and is expected to become large at the steepest edges, hence the term 
"edgeness" . The quantities u a and nothing but the trace and determinant of 

the Hessian matrix of the field's profile and thus provide information about the nature of extrema. 
Finally, related to the curvature of isolines, which becomes large at corners 



of the field's isolines and was accordingly named "cornerness" by ter Haar Romeny et al. (1991). 

In Figure we demonstrate the effect of scale-space operators and geometric invariants on a simple 
pattern. 

1 We use the summation convention for pairwise indices in a product, and denote partial derivatives by indices 
following a comma. 
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2.3 Statistical Descriptors 



After applying scale-space filters and calculating geometric invariants, we arrive at a field X a that 
contains less unwanted structure for an appropriate choice of a and enhances geometric features 
according to the invariant that has been chosen. Several possibilities exist to construct statistical 
descriptors from these fields. 



2.3.1 Isodensity contours 



The surface integration necessary for the calculation of the genus of isodensity contours flGott III 



et al. 1986| ) can be reduced to taking a spatial mean value of Koenderink invariants, multiplied 



with a delta function for selecting the specific isodensity contour. For the Euler characteristic 
X = 1 — g of the isodensity contour of the field u{x) to the threshold v, for example, one obtains 

27T J UhUh 



Recently, this rationale has been used by (Schmalzing fc Buchert (1997) to calculate not only the 



genus, but all Minkowski functionals for a continuous random field in three dimensions. 



2.3.2 Excursion sets 

It is possible to measure the size of an invariant's excursion set over its mean value - the sharper 
the features emphasized by the invariant, the smaller the excursion set. Double Poisson processes 
for various geometrical objects can be efficiently discriminated by measuring the excursion sets 
and successively filtering out larger and larger scales in the point distribution (see |Schmalzing| 
1996| for further details). 



2.3.3 Probability density 

Finally, one can think of using the probability density of an invariant as an indicator of certain 
patterns. It is especially promising to compare the measured curves to analytically calculated 
mean values for simple types of random fields. We have developed a test for the Gaussianity of 
random fields and will present preliminary results of its application to CMB anisotropics in the 
following Section 9. 



3 Gaussianity of Random Fields 

It is generally believed that the Cosmic Microwave Background anisotropics obey the statistics of 
a Gaussian random field, at least on scales observed by the COBE satellite. On smaller scales the 
situation is not so clear. In any case it is important to assess whether there are any non-Gaussian 
signatures in the microwave sky. 



3.1 Gaussian Random Fields 



The statistical properties of homogeneous and isotropic Gaussian random fields are solely de- 
termined by their power spectrum P(k) or equivalently their two-point correlation function £(r). 
Among other works, the standard textbook by Adler (1981) or the article by Bardeen et al. (1986)| 
are devoted to an extensive study of random fields in general and Gaussian ones in particular; the 
calculational schemes presented there also allow to tackle Koenderink invariants of random fields 
analytically. 
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Figure 3: This figure shows a Microwave sky map constructed from the COBE DMR four-year 
data. By removing the monopole and dipole from the mean value (A-\-B)/2 of the 53GHZ channels, 
the anisotropics on smaller scales become visible. Finally, we only use the 3881 pixels not covered 
by COBE group's mask for galactic emission to plot this map and to apply Koenderink invariants. 



3.2 Koenderink Invariants 



It turns out that the probability densities of any Koenderink invariant (see Schmalzing 1996 for 
detailed results) depends solely on three common parameters (£(0), £"(0) and ^""(O)). Further- 
more, they only determine the width of the distribution, while the shape is the same for any 
Gaussian random field. They can be fitted using the functions measured from the random field, 
and the goodness of the fit gives a measure of the deviation from Gaussianity. 

The panels shown in Figure || illustrate the method outlined above. We show a random field, 
the CO BE (A+B)/2 signal in the 53GHz channel with a customized galactic cut (details can be 



found in [Bennett et al. 1996[) . From this random field several distribution functions of Koenderink 
invariants have been calculated. Although the behaviour is largely consistent with the assumption 
of a Gaussian random field, some deviations are visible. Figure || shows and explains one of the 
seven panels in more detail. 



4 Outlook 



Apart from the example discussed in this talk, Koenderink invariants suggest a variety of further 
studies in cosmology. 



4.1 Large Scale Structure 

The possibility of extending Koenderink invariants to arbitrary dimensions has received only brief 
attention in this article. On the one hand, Koenderink filters can be used to construct morpho- 



logical statistics for the three-dimensional matter distribution in the Universe ( Schmalzing 1996 ), 
on the other hand, the insights gained from the theoretical fou ndations may be used to shed new 
light on seemingly unrelated measures such as genus statistics (|Schmalzing fc Buchert 1997 ). 



4.2 Microwave Background 

The application discussed in this talk uses Koenderink invariants to assess whether non-Gaussian 
features can be seen in the Cosmic Microwave Background. The preliminary results suggest that 
this is not the case, at least for the scale probed by the COBE satellite. However, a number of 
issues remain to be addressed. 



Numerical and, if possible, analytical calculations for certain non-Gaussian fields can test the 



Probability Density 




Figure 4: Here we display the results of an analysis of the map in the previous figure. The 
probability densities of several Koenderink invariants are shown, both as measured from the data 
(solid line) and as calculated for 100 realizations of a Gaussian random field (the shaded area 
denotes la fluctuations). 
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Probability Density 








Figure 5: This plot shows the probability density of one of the invariants from Figure |. The solid, 
irregular line was measured from the COBE DMR data. The area indicated by the three smoother 
solid lines corresponds to la-fluctuations calculated from one hundred realizations of a Gaussian 
random field. Obviously the fluctuations of the data around the fitted mean value (dotted line) 
are consistent with a Gaussian random field. However, we do see slight discrepancies; future work 
will assess whether they become significant at higher resolution, when improved statistics reduces 
the sample-to-sample variance. 
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significance of our findings and test the discriminative power in comparison to other methods (e.g. 
Kogut et al. 1995| ). 



With regards to the forthcoming high-resolution surveys of the microwave sky ( Bersanelli et al 
1996| ; Bennett et al. 1995| ), the performance of Koenderink filters on smaller scales needs to be 



tested. Their noise reduction abilities are another important subject of study. 

Finally, it would be interesting to apply related descripto rs o ther than the probability densities of 
Koenderink invariants. The methods outlined in Section |2.3| are promising candidates. 
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